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^ . Abstract 

We apply an Osterwalder-Seiler version of twisted mass QCD to a study of the Bx pa- 

" ^ I rameter, in which three of the four quark fields making up the relevant AS" = 2 operator 

' are maximally twisted with the same twist angle, while the fourth one has a twist angle 

■ of opposite sign. It is known that this setup ensures automatic improvement of the bare 
, K^-K^ operator matrix element and multiplicative renormalization of the AS" = 2 oper- 
\ ator, at the price of breaking the K^-K mass degeneracy by discretization effects. As a 
' result, two dominant systematic errors of the -Bk determination with Wilson fermions are 

kept under control. With the Clover term included in the fcrmion action, we perform a 
. . , feasibility study and find, in the quenched approximation , a significant improvement of the 

^ ' scaling behaviour of i?Kj compared to earlier standard tmQCD determinations. Moreover, 

' we study in detail the K^-K mass splitting that characterizes this approach and confirm 

^ . that, in the presence of the Clover term, it is greatly reduced in a maximally twisted theory. 



^On leave of absence from DESY, Platanenallee 6, D-15738 Zeuthen, Germany. 



1 Introduction 



The bag parameter of neutral iC-meson oscillations B^, has been the object of 
many lattice QCD computations. Several discretizations of lattice fermions have 
been implemented for this purpose. For the most recent quenched result see ref. [1]. 
Attention has now shifted to unquenched estimates; see ref. [2] for a recent review. 

The Wilson fermion results of Bk are normally the least accurate, due to a 
limited control of two sources of systematic error: 

(i) Loss of chiral symmetry causes the relevant AS" = 2 four-fermion operator to 
mix, under renormalization, with four other operators of the same dimension 
and different chiral representation [3]. In lattice discretizations which do not 
violate chiral symmetry this operator is multiplicatively renormalizable. 

(ii) Discretization effects of the Bk estimate at finite lattice spacing are 0{a); 
with staggered and Ginsparg- Wilson fermions they are O(a^). The traditional 
remedy of Symanzik improvement is not viable here, because of the large 
number of higher dimensional 0{a) countcrtcrms required. 

The first of these problems has been dealt with in ref. [4], by using Ward 
identities, and in refs. [5,6], by implementing twisted Wilson fermions [7]. In the 
latter case, the lattice fermion action consists of a twisted up-down fermion doublet 
(tmQCD) and a standard (untwisted) Wilson strange fermion. A second possibility, 
valid only in the quenched approximation, is that of a tmQCD discretization of 
a degenerate down-strange doublet. The lack of mixing with other operators in 
both tmQCD variants is achieved through the mapping of the usual, parity-even 
AS = 2 operator to its parity-odd counterpart, which is known to be multiplicatively 
renormalizable, even with Wilson fermions [8,9]. Both formalisms have been applied 
in refs. [5,6], combined with the Schrodinger functional renormalization and Re- 
running of ref. [10] . 

In order to implement 0(a) improvement, without the use of Symanzik coun- 
terterms, one must ensure that all flavours are regularized in a tmQCD framework at 
maximal twist. In other words, the twist angle for all quark flavours must be tuned 
to ±7r/2, while the mapping from parity-even to parity-odd AS = 2 operator must 
be preserved. As discussed in ref. [11] this is not possible in the standard tmQCD 
formalism. A way round this problem was provided in ref. [12], by using somewhat 
different regularizations for sea and valence quarks. The former may be standard 
tmQCD, while the latter is the so-called Osterwalder-Seiler variant of tmQCD. Each 
valence quark is maximally twisted to ±7r/2. However, unlike in standard tmQCD, 
the valence flavours are not combined into isospin doublets. The four-fermion op- 
erator consists of four distinct, maximally twisted valence flavours, three of which 
have twist angle with the same sign (say, +Tr/2) while the fourth twist angle has 
the opposite sign (say, — 7r/2). This combination of valence quarks ensures that 
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the AS = 2 operator of interest is multiplicatively renormalizable and its — K 
matrix element is automatically improved; the same is true for 5k. 

It is not surprising that a price has to be paid for these advantages. In the 
setup described above, the i^'^-meson consists of a strange/down valence quark- 
antiquark pair with, say, the same twist angle, while the K -meson has a pair with 
opposite twist angles. This means that the two mesons have (pseudoscalar) masses 
which differ by O(a^) discretization effects. Although in principle this mass splitting 
vanishes in the continuum limit, it may be quite sizeable at finite lattice spacing. 
In the quenched approximation such effects have indeed been studied in the past 
(cf. refs [13-15]), being indirect manifestations of flavour symmetry breaking by 
tmQCD. 

In this work we have also performed a similar study of such flavour symmetry 
breaking effects for the ivT-meson mass mx and decay constant /k. A comparison of 
our results with those of refs. [13, 14] indicates that the presence of the Clover term 
in the action greatly reduces such flavour breaking systematics, in accordance with 
the findings of an earlier study [15]. Moreover, we find that the scaling of i?K with 
the lattice spacing is significantly improved, if compared to the previous standard 
tmQCD quenched studies of refs. [5,6]. 

2 Twisted and Osterwalder-Seiler valence quarks at maximal twist 

The setTip of our formalism follows very closely that of ref. [12]. Each valence quark 
field qf (the subscript labels flavour, = u,d,s,...) is discretized by a fermionic 
action of the form 

Sf = a'^Y^ qf{x)[D^ + mf +ij5iJ,f]qf{x) , (2.1) 

X 

where is the standard Wilson-Dirac fermion matrix with a Clover term. The 
Wilson plaquette action is the regularization of the pure gauge sector of the the- 
ory. In this work all flavours are degenerate; this is only a matter of choice. As 
our simulations are performed in the quenched approximation, wc will ignore is- 
sues concerning the lattice regularization of the sea quarks; see [12] for a related 
discussion. 

Let us consider quark bilinear operators with distinct flavours. We distinguish 

two cases, depending on the relative sign of the twisted mass terms (or equivalently, 
the twist angle) of these flavours. For example, the pseudoscalar density 

PdTix) = d{x) 75 s{x) (with fid = -M.) (2.2) 

is said to be of "twisted mass" (tm) type, if the corresponding twisted mass terms 
have opposite signs (i.e. fij, = —Us)- In tbis case the two flavours may be grouped 
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in a doublet, with a mass term of the form ifJ^^^r^, corresponding to the standard 
tmQCD formahsm (rs is the Pauh isospin matrix). A second possibihty is that 
of lid = Ms; this will be referred to as the Osterwalder-Seiler (OS) case [16]. The 
pseudoscalar density is then denoted by 

PrfOS(x) = d»75s(x) (with = ;u,) . (2.3) 

Analogous definitions hold for other quark bilinear operators. The OS discretization 
is only applicable to valence quarks. Sea quarks, if regularized by the tmQCD lattice 
action, must be organized in doublets, with a flavour non-singlet twisted mass term 
i//75r3, in order to avoid the generation of an unwanted O-teim. The implementation 
of OS valence flavours is a mixed action formulation. Inevitably, in full (unquenched) 
QCD this introduces unitarity violation at finite lattice spacing. 

The specific continuum operators of interest to us are Pds{x) = d{x)^5s{x) and 
-Afi^si^) — d{x)"f ^'^^s{x) . As shown in ref. [7], once the lattice theory is renormalized 
in a mass-independent renormalization scheme, such quantities are related, up to 
cutoff effects, to those of the lattice twisted theory through the quark field chiral 
rotation^ 

i^'°^\x) = exp[l75r3a] [V(x)]r V^'°°*(x) = ^{x)]^ expffTsTsa] ; (2.4) 

the phase a is the so-called twist-angle, defined through tan(a) = /XR/mR. This of 
course refers to the standard tmQCD formulation, with an isospin quark doublet ijj 
and a twisted mass term of the form inipj^T^ip. For the case under consideration, 
with a distinct action (|2.ip for each fiavour qj, the corresponding rotations of the 
valence quark fields are 

qjixr^' = exp[i75a^] [g/(x)]R g/(x)™'^* = [g/(x)]R exp[i75a^] , (2.5) 

where the sign of the phase «/ is that of the corresponding mass term /xj . In order to 
ensure automatic improvement, we consider the case of maximal twist (a/ = ±7r/2), 
for which the above chiral rotations induce the following relations for the two-fermion 
operators under consideration: 

(2.6) 
(2.7) 
(2.8) 
(2.9) 

All such relations, understood to be shorthand expressions for equations between 
correlation functions (or matrix elements) involving these operators, are valid up to 
discretization effects. 



pcont 
^ds 




= ZpP'^ 


/I cont 




= -iZvV^4s 


pcont 
^ds 




= iZsSZ"" 


^cont 
■^u,ds 


= [^?i]R 


7 /lOS 



^Our notation is the following: the superscript "cont" denotes continuum quantities, while [■ 
stands for lattice (re) normalized quantities, corresponding to either tm or OS discretization. 
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Maximal twist is achieved by tuning the standard mass parameter mo to its 
critical value rricr (with corresponding value Kcr of the hopping parameter). There 
are several ways of doing this. A popular procedure consists in choosing a two-point 
correlation function which in the continuum violates parity and flavour symmetry 
and must therefore vanish. Thus, at fixed fi, the point mo = mcr(M) is determined, 
for which the bare correlation function vanishes. Finally, the critical mass rricr is 
found by extrapolating mcr(/i) to = (unless the smallncss of /i renders the 
extrapolation unnecessary). For details see ref. [17]. This procedure, based entirely 
on a tmQCD setup, ensures the absence of 0{a?'^) discretization effects (with k 
integer) in the determination of rricr, while 0{a) effects are present. On the other 
hand, physical quantities like hadronic masses and matrix elements are automatically 
0(o)-improved; more precisely all 0{a'^^^^) discretization effects are absent. 

In the present work, we follow a very different approach, described in detail 
in refs. [5,6]. It consists in using the standard Schrodinger functional estimate of 
rricr, by requiring the vanishing of the PCAC quark mass, as in ref. [18]. This pro- 
cedure uses standard (non-twisted) Wilson fermions, thus obtaining an explicitly 
/x-independent estimate of rricr- However, in order for this rricr to be 0(a)-improved, 
both the fermionic action and the axial current in the PCAC relation must have 
0{a) Symanzik-counterterms. Therefore the non-perturbative estimates of csw and 
CA, determined with standard Wilson fermions are essential in this procedure. With 
rricr thus obtained in a non-twisted, improved Wilson fermion setup, we compute 
hadronic masses and matrix elements in a maximally twisted framework. This ap- 
proach has also been adopted in ref. [15], with lattices having periodic, rather than 
Schrodinger functional, boundary conditions. In that work, a detailed comparison 
is made, at /3 = 6.0, between pseudoscalar masses and decay constants (of the tm 
variety), obtained with the two estimates of Kcr described above. 

An important consequence, shared by these methods of tuning to maximal twist, 
is that the so-called "chirally enhanced" discretization effects of the form a^/m^, 
present in lattice correlation functions of fermion multi-local operators, are reduced 
to X /m^\ cf. ref. [17] . This is however of little relevance to our results obtained 
close to the iC-meson mass region, i.e. far from the chiral limit. 

3 Results 

We performed quenched simulations at three values of the gauge coupling, namely 
(3 = 6.0, 6.1, 6.2. At each coupling we have tuned the degenerate quark masses 
a/j.^ = ajJLs to a couple of values which are close to half the strange quark mass, so 
as to simulate charged K-mesons with degenerate flavours. In practice this entails 
tuning the twisted mass values so as to obtain pseudoscalar meson masses close to 
that of the physical Kaon. The values of Kcr, used in the simulations in order to 
achieve maximal twist, are collected in Table 3 of ref. [6]. We use the renormalization 
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constants, improvement coefficients and the ratio r^ja as collected in Appendix A 
of ref. [5]. Physical quantities are usually expressed in units of tq in the present 
work [19]. 

The details of our simulations are listed in Table [TJ The lattice size is L'^ x T 
with standard Schrodinger functional boundaries. The lattice calibration, expressed 
in terms of a/2ro, is that of ref. [5]. In the Table we also show the range (in units of 
xqIItq) for which the ground state of our correlation functions (i.e. the K-meson) 
has been isolated. 





a 
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L 2rn ' 2rn J 
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(24, 
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[1.30, 3.17] 
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0.0135 
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0.0931 


(24, 
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[1.30, 3.17] 
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(24, 
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[1.26, 3.47] 


100 


0.0125 


I 


6.1 


0.0789 


(24, 


60) 


[1.26, 3.47] 


100 


0.0110 


II 


6.2 


0.0677 


(32, 


72) 


[1.35, 3.52] 


50 


0.0105 


I 


6.2 


0.0677 


(32, 


72) 


[1.35, 3.52] 


50 


0.0090 


II 



Table 1: Details of the runs for three values of the gauge coupling. 



3.1 Comparison of discretization effects between tm and OS quantities 

All observables computed in this work are obtained from the large time asymptotic 
limit of operator correlation functions with Schrodinger functional boundary con- 
ditions. The notation is standard, following closely that adopted e.g. in ref. [5]. 
For instance, /a (/a) denotes the Schrodinger functional correlation function be- 
tween a fermionic operator Aq in the bulk and a pseudoscalar boundary operator 
CslbCd at time-slice = {xq = T). All such correlation functions are properly 
(anti)symmetrized in time, when used to extract effective pseudoscalar masses and 
decay constants. These quantities (as well as in subsection 13. 2p are obtained at 
large time-separations from the boundaries, in order to avoid contamination from 
higher excited states. Moreover, in this asymptotic limit, 0{a) discretization effects 
from the boundaries may be considered negligible. 

We have measured pseudoscalar effective masses aM(xo), both in the tm and 
OS framework: 



aMt-(xo) = ^ln 



f}f^{xo + a) 



aM°S(xo) = ^In 



(3.1) 
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1.169(5) 
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17 


6.1 
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1.295(6) 


1.349(8) 


8 


6.1 


II 


1.222(6) 


1.279(8) 


9 


6.2 


I 


1.268(8) 


1.301(10) 


5 


6.2 


II 


1.182(8) 


1.217(10) 
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Table 2: Pseudoscalar masses of tm and OS type and their mass splitting as defined 
in Eq 



In Table [2] we show the results of the pseudoscalar meson masses from both reg- 
ularisations. Recall that the physical value of the Kaon mass in units of tq is 
to-Mk = 1.2544. Therefore, for all three /3-values, the K-meson decay constant Fk 
and the bag parameter i?K discussed below, are obtained through short "interpola- 
tions" between two points. We also show in Table [2] that the relative mass splitting, 
defined by 

Rm ^:^M^2 ' (3-2) 

decreases with decreasing lattice spacing. This is a clear evidence that the mass 
splitting is a discretisation effect. 

In order to compare our results to those of other collaborations [13-15], we 
consider the pseudoscalar meson mass splitting in units of rg: 

A(r^M2) = [roM^^f - [r^M^'^f . (3.3) 

The tm-OS pseudoscalar mass splitting may be affected by two factors: (i) the 
presence of the Clover term and (ii) the way Kcr is determined in order to ensure 
maximal twist. We first recall the different choices of fermion action and Kcr, made 
in refs. [13-15]: 

• Both refs. [13, 14] have data obtained with a maximally twisted Wilson ac- 
tion [without a Clover term) but implement different methods for computing 
Kcr- In particular, ref. [13] uses two methods, namely the PC AC relation for 
maximally twisted quarks (yielding an "optimal" k,ct) and the vanishing of the 
pseudoscalar meson mass in the standard (untwisted) theory. On the other 
hand the authors of ref. [14] choose to tune Kct through the vanishing of the 
parity-odd correlation function (y^°"*'(a;)i-*™'^*(0)), corresponding to another 
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Figure 1: Square- mass splitting between tm and OS pseudoscalar mesons, as a 
function of the tm pseudoscalar meson mass: (a) Comparison between our results 
and those of ref. [13]; the latter are obtained using the PC AC determination of Kqt- 
(b) Comparison between our results and those of ref. [14]. 



definition of an "optimal" A further difference is that in ref. [14], the values 
of Kcr(/^), computed at each twisted mass, are not extrapolated to /i — > 0. 

• In ref. [15] the Clover term is included in the fermion action, just like in the 
present work. Their determination of Kcr is based on the PCAC relation, eval- 
uated in the standard Wilson theory (i.e. zero twist) with periodic boundary 
conditions, at large pseudoscalar masses, extrapolated to the chiral limit. On 
the other hand, our Kct determination is based on the PCAC relation with 
Schrodinger Funtional boundary conditions, extrapolated to the chiral limit 
from small masse^. Both estimates are 0(a)-improved. 

Comparisons between results of these works are mostly limited to pseudoscalar 
mesons in the region of Mk- Pseudoscalar masses in the lighter mass range may 
be subject to finite volume effects (e.g. L = 16^ at /? = 6.0), while those related to 
heavier masses are subject to larger discretization errors. A close look at the details 
of these works suggests that the tm pseudoscalar masses in the Mx range agree 
reasonably well, while the OS ones do not. Thus any discrepancies in A(rQM^) are 
due to the latter mass. We may therefore conclude that these works indicate that 

^In our simulations, the twisted mass angle, measured from the ratio dofv^/dofA^, varies be- 
tween 87° -90°. 
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Figure 2: Square-mass splitting between tm and OS i^-mesons, as a function of tlie 
lattice spacing squared. The data of ref. [13], displayed here, are obtained using the 
PC AC determination of Knr. 



the pseudoscalar mass M*™ is weakly dependent on the presence or absence of a 
Clover term and/or the way maximally twist is implemented. 

Next, we examine the influence of Kcr in the evaluation of M^^ and consequently 
of A(rQM^). Our conclusions may be summarized as follows: 

(i) The results of ref. [13], obtained with two choices of Kcr, indicate that, for the 
coarser lattices, these quantities are affected by the way Kcr is determined. 

(ii) The comparison made in ref. [14], between their M^^ and those of ref. [13], 
obtained with the PC AC determination of Kcr, also shows discrepancies (cf. 
Fig. 1 of ref. [14]; see also our Figs. [H [2]). Ag ain the choice of affects AfOS 
and A(rgM2). 

(iii) We observe a better agreement between A(rgM^) of ref. [14] and that of 
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ref. [13], obtained with the pseudoscalar determination of Kcr- 

(iv) From Fig. 5 of ref. [15], we estimate that at the ii'-meson region, A(rgM^) 
0.27 (for P = 6.0), while in our computations we find A(rgM^) ~ 0.25. 

The good agreement of the latter estimates suggests that once the Clover term is 
included in the action and an 0(a)-improved Kct is used to tune to maximal twist, 
the tm-OS pseudoscalar mass splitting is only mildly affected by the details of the 
Kcr determination. Conversely, it appears that the 0{a) effects in the "optimal" 
determination of Kcr (with a Wilson action without a Clover term) may somehow 
induce large O(a^) effects in M^^ and consequently in A(rQM^). 

We now investigate the direct influence of the Clover term in the mass splitting 
A(rQM^). In Fig. [1] we plot our data against those of refs. [13, 14]. The mass 
splitting of our results is sig nificantly smallefl We attribute tkis to the inclusion 
of the Clover term in the action. This effect, already seen in ref. [15] at /5 = 6.0, is 
confirmed here for more lattice spacings. As the data of refs. [13,14] clearly indicate, 
this splitting depends weakly upon the quark mass values. We expect a similarly 
weak mass dependence in our case. 

In Fig. [2] we plot A(rgM^) against the lattice spacing squared, for the present 
work and for refs. [13,14]. In all cases this cutoff effect is decreasing with the lattice 
spacing. Contrary to expectations, there is no clear evidence for its vanishing in 
the continuum limil[f|. This is probably due to the fact that we only have results at 
three lattice spacings for each case. Moreover, our lattices are probably too coarse 
to reveal the vanishing of A(rQM^) in the continuum limit|f|. 

Finally, we point out that this mass splitting has been recently studied in an 
unquenched Nf = 2 framework, with a maximally tmQCD action without a Clover 
term; see ref. [20]. In this preliminary work, pronounced Rm values of up to about 
50% have been reported. 

A similar analysis is performed for the pseudoscalar meson decay constants. 
In the Schrodinger functional framework they are obtained from the axial current 
correlation functions /aj^, properly normalized by the boundary-to-boundary corre- 
lation function /i; see [21] for details. For the maximally twisted tm and OS cases 

^In refs. [13, 14] the pseudoscalar masses aA/'™ and aM'~'^ are tabulated. From these we can 
easily compute A(roM^). The error on the latter quantity has been read off from Fig. 2 of ref. [13] 
and Fig. 3 of ref. [14]. This rough error estimate is adequate for the present qualitative comparison 
with our data. 

^This contradicts the claim of ref. [13], where the mass splitting is shown to vanish in the 
continuum, for lighter pseudoscalar massess. In this respect we wish to point out that the tm and OS 
pseudoscalars were computed in that work with much bigger errors. Moreover, the L = 16^, /3 = 6.0 
results at the lightest values of afi may not be free of finite size effects. Thus we consider the 
conclusion of that work on the mass splitting under consideration as not definitive. 

®Most peculiarly, the vanishing of A(roAf^) in the continuum appears to be supported by our 
data if plotted against a**. The data is well fit either by a function of the form Co + CiO^ or by one 
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under investigation, the specific expressions for the decay constants F*™ and F^^ 
are (in the large-time asymptotic regime): 

F*'" ^ 2(M^'^L=^)-^/^ exp[(xo - T/2)M^'^] , (3.4) 

^ 2(Af QSl^)-V2 exp[(xo - r/2)MQS] ^^/a^o) . (3.5) 

In practice the above quantities are obtained in the xo-range in which the pseu- 
doscalar effective masses have been extracted. 

A further method for computing F is based on the PCVC relation, expressed 
in terms of Schrodinger functional correlation functions: 

- M*-Zy/r (xo) = 2i/i/*'"(xo) . (3.6) 
The corresponding decay constant is computed as 

fPCVC ^ -4^(M*-L3)-V2exp[(xo-T/2)M*-]^^^. (3.7) 

In Table [3] we show our results for the pseudoscalar decay constants. These 
values are in a general agreement with the corresponding estimates in Tables 6 and 
7 of ref. [6], but the errors are bigger, due to reduced statistic^. The principal 
observation is that F*™ and F^^ are compatible within errors. 



3.2 Improved i?K parameter 

We now pass to the computation of the i?K parameter, on the lines of ref. [12], which 
ensures that the bare i?K is automatically improved. The four-fermion operator of 
interest is 

Qvv+AA = [sif^d] [sj^^d] + [s-f^j5d] [57^750!] ■ (3.8) 

On the lattice, we discretize this operator as proposed in ref. [12]; i.e. it is expressed 
in terms of current products made of both tm and OS type, at maximal twist: 

QTl+Iv = [si.dri^^.l^df'' + [si,Todr[s^,df^ . (3.9) 

Using the chiral rotations of eq. (j2.5p . we obtain the relation between the lattice and 
continuum operators, up to discretization effects: 

/Ocont rA-)tm— OS i ry /^tm— OS /o i n^ 

'-<iVV+AA — [VvA+AvJR — ■^VA+AV'^VA+AV [6.1U) 

''The quantities roF^''"^^ and ro-F*™ have also been computed in ref. [6]. However these com- 
putations had been performed at significantly heavier quark masses and extrapolated to the Kaon 
mass value. This accounts for the 4% difference between these results and ours at /3 = 6.0; this 
difference decreases with increasing /3. 
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Table 3: Pseudoscalar decay constants obtained from PCVC, and the axial current 
of tm and OS type. We also show the interpolations at the physical K-meson point. 



The main reason behind using this mixed tm-OS formalism is the fact that in this 
way the -Bk matrix element, computed with (twisted) Wilson fermions, is both 
automatically improved and multiplicatively renormalizable. 
-Bk is obtained form the ratio 

■ ry 77tm — OS 

where -^^'y^+^y is the Schrodinger functional correlation function with the four- 
fermion operator of eq. (|3.9p in the bulk and the usual boundary sources at the time 
edges. Note that these sources are of tm type at time-slice xq = and of OS type 
at time-slice xq = T. Thus, away from the time boundaries, the correlation function 
F^'X+Jiv displays an asymptotic behaviour of the form exp[— M*™xo] exp[— M*-'^(r — 
a^o)]. In order to match this choice and cancel the exponentials, the denominator 
consists of a correlation function fy^ of the tm type, involving the xq = boundary, 
and a correlation function of the OS-type, involving the xq = T boundary. The 
axial and vector currents in these correlation functions are the corresponding tm 
and OS twisted versions of the axial current; cf. eqs. (j2.7p . (j2.9p . In our analysis 
we have also considered the symmetric situation with tm and OS sources exchanged 
and the denominator of the above equation adjusted accordingly. The -Bk result 
quoted below is the average of the two estimates thus obtained. The quality of our 
raw results for Rbk is shown in Fig. [3l 

At large time separations from the boundary, the bare i?K estimate is automati- 
cally improved. The renormalization constant ZyA+AV is known non-perturbatively 
from ref. [10]. As it is computed in the chiral limit but is not improved, it suffers 



11 



1.6 



0.8 



0.6 




^ \ ' \ ' — 

■ P=6.0, dataset I (shifted by +0.2) 
A |3=6.1, dataset I (shifted by +0.1) 
• |3=6.2, dataset I 



0.2 



0.4 



0.6 



0.8 



Figure 3: The ratio R^hgi as a function of time, for three values of the gauge 
coupling. The data for /3 = 6.0 and 6.1 have been shifted for clarity. The horizontal 
lines indicate the i?K averages and the extension of the plateau. 



from O(aAQCD) discretization effects. These are not expected to dominate the scal- 
ing behaviour of Sk! evidence for this will be provided below. 

In Table H] we list our B^^^ estimates. The scaling of these results is to be 
compared to the ones of refs. [5,6]. We remind the reader that the latter had been 
computed in two distinct frameworks: (i) the so-called 7r/2 case, consisting of a 
maximally twisted down quark and a non-twisted strange one and (ii) the so-called 
7r/4 case, in which down and strange flavours were in the same isospin doublet with 
a twist angle 7r/4. Both frameworks have quark flavours which are not maximally 
twisted and this implies that the results of refs. [5,6] are not improved. 

In this respect, we wish to draw the reader's attention to a subtlety, related to 
the Symanzik-improvement term acAd^P of the axial currents in the -Bk denomi- 
nator. In these earlier works cited above, the influence of this term in the scaling 
behaviour of i?K was carefully monitored. The result of this analysis is summarized 
in Fig. IHa). Note that scales much better when = (i.e. when the i?K 

denominator is not improved), than when it is tuned so as to eliminate 0{a) effects 
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(3 dataset B] 



6.0 


I 


0.772(10)(9) 


D.U 


TT 
ii 


(J. (b6{lZ)\\3) 






(J. / /Z(iL)j(yj(i4j 


0.1 


T 
i 


n 7(^0/' 1 0^ /'o^ 

U. (KyAylA )\\)) 


6.1 


II 


0.746(13)(9) 




toMk 


0.753(12)(9)(15) 


6.2 


I 


0.758(12)(9) 


6.2 


II 


0.742(13)(9) 




toMk 


0.756(12)(9)(15) 



Table 4: results of the present work. The errors are, in order of appearance: 

(i) statistical, (ii) uncertainty due to the renormalization constants ZyA+AV-^ 
and Zy^ (iii) total error, obtained by adding the first two in quadrature. 
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Figure 4: (a) Scaling behaviour of the results of refs. [5,6] for: (i) the 7r/2 case 
with a CA term (empty squares) and with = (filled squares); (ii) the 7r/4 case 
with a CA term (empty diamonds) and with Cyi = (filled diamonds), (b) Scaling 
behaviour of the results from the present work with Cyi = (filled squares) 

and with a ca term (empty squares). Symbols have been shifted for clarity. See the 
text for details on the continuum limit extrapolations (C.L.). 



from the axial current. In ref. [5] it was speculated that this could be due to some 
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cancellation between the discretization effects of the -Bk numerator and denomi- 
nator, which is spoilt once the ca_ term is switched on. Whatever the reason, the 
data of Fig. IHa) clearly indicate that 0{a) effects influence significantly the scaling 
behaviour of Sk- The continuum limit value shown in this figure is the final result 
quoted in ref. [6]. It has been obtained by a combined fit, linear in a/2ro, of the 7r/2 
and 7r/4 data with a ca_ term in the denominator (which is the most conservative 
option) . 

In Fig. m^b) we show our results, plotted against the linear lattice spacing a/2ro. 
They are shown both for ca = and for non-zero ca, in order to monitor the influence 
of the relevant counterterm. Note that due to automatic improvement, this is now an 
O(a^) effect. The compatibility between the two sets of results in Fig.|Hb) indicates 
that, unlike the previous case, the ca term does not spoil the good scaling properties 
of Bk- Moreover, the almost fiat scaling behaviour may suggest that the O(aAQCD) 
discretization effects of Zva+AV: as well as the 0(/i^a^) effects of the — 
matrix element, are small. In the light of this we have performed two continuum 
limit extrapolations, one linear in a and another linear in a^. Both extrapolation 
results are shown in Fig. [Hb)), for the ca = data. 

The linear fit of our B^'^^ data with respect to a gives the following continuum 
result: 

= 0.706(65) ix^/dof = 0.30) 
while the linear fit with respect to gives: 

S^Gi ^ 0.733(34) ix^/dof = 0.26) . 

Both results compare nicely with the published value: B^'^^ = 0.735(71) of ref. [6]. 
Our error in the latter value is smaller because the extrapolation in is shorter. 

4 Conclusions 

In this work we have addressed a couple of issues related to the proposal of ref. [12] for 
the computation of i?K in a tmQCD setup, using a four-fermion operator made of two 
tm-type quarks and two OS-type quarks, all maximally twisted. This ensures that 
the operator is multiplicatively renormalizable and automatically improved. The 
price to pay is that in this way the Kaon and anti-Kaon states of the weak matrix 
element {K \0^^^'^\K'^) are composed of different lattice quark field combinations 
(tm and OS) and thus are non degenerate at finite lattice spacing. 

We have performed a quenched study in order to address these issues. In par- 
ticular, we find that the inclusion of the Clover term in the fermionic action sig- 
nificantly reduces the mass splitting between Kaon and anti-Kaon at finite lattice 
spacing. This result, already observed in an earlier work at a single lattice spacing, 
has been confirmed by us for several /3-values. It should also be noted that in earlier 
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studies, in which a Clover term was not included in the action, this mass splitting, 
besides being much bigger than the one we measure, was found to be significantly 
affected by the choice of non-perturbative procedure for tuning the lattice theory 
to maximal twist. This does not appear to be the case once an 0{a) Symanzik- 
improved procedure, with a standard (untwisted) mass term, is used to tune the 
hopping parameter to its critical value. We also confirm that this splitting decreases 
as the continuum limit is approached. 

The second objective of the present paper is the presentation of first results for 
Bk, obtained in this tm-OS setup. We have compared our results to the earlier ones 
of refs. [5,6], in which i?K had been computed in a tmQCD setup, without maximal 
twist on all flavours and thus without automatic 0(a)-improvement. These results 
were rather sensitive to discretization errors, as revealed by the dependence of Sk 
on the Symanzik 0{a) term of the axial current at coarser lattice spacings. In the 
present work, bare matrix elements are automatically 0(a)-improved, leaving us 
only with O(aAQCD) effects from the renormalization constant of the four-fermion 
operator. Our results display significantly better scaling, compatible with an O(a^) 
behaviour. 
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